NUCLEAR SCIENCE AND TECHNIQUES 26, 050402 (2015) 


New methods to remove baseline drift in trapezoidal pulse shaping* 


HONG Xu (#E}IH),! MA Ying-Jie (3375), t ZHOU Jian-Bin (H 2250),! CHEN Tie-Guang (fRFKIL),” 
HU Yun-Chuan (#54 ])!|),! WAN Xin-Feng (77 #1Il),! and DU Xin GÆ)! 
‘College of Nuclear Technology and Automation Engineering, 
Chengdu University of Technology, Chengdu 610059, China 
°Chengdu Gaotong Isotope CO., LTD. (CNNC), Chengdu 610000, China 
(Received May 8, 2014; accepted in revised form August 18, 2014; published online October 20, 2015) 


Trapezoidal pulse shaping algorithm is widely applied to improve signal-to-noise ratio (SNR), throughput 
and energy resolution with the properties of noise suppression, pile-up pulse separation and ballistic deficit 
correction. The algorithm can be acquired by z transform method which is easier for derivation. However, the 
baseline drift of trapezoidal pulse appears because the noise superimposes on the input signal. In this paper, 
two new methods based on convergence analysis and noise suppression are proposed to remove the baseline 
drift resulting from trapezoidal pulse shaping. Simulations and experimental tests are carried out to verify the 
methods. The results demonstrate that the proposed methods can remove baseline drift in trapezoidal pulse 


shaping. 
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I. INTRODUCTION 


Digital pulse shaping algorithm is a main contents of dig- 
ital spectrometer. They are injected to realize pulse shaping 
in field programmable gate array (FPGA). Trapezoidal pulse, 
with its rise time being equal to the fall time, provides a near 
optimum signal-to-noise ratio (SNR). The flat top of trape- 
zoidal pulse can be flexibly set to adapt different measure- 
ment conditions. Therefore, trapezoidal pulse shaping algo- 
rithm is widely applied to improve SNR, throughout and en- 
ergy resolution. 

The theory of trapezoidal pulse shaping has been well de- 
veloped and published. Radeka designed a trapezoidal filter 
based on a gated filter system, which achieved good resolu- 
tion with large germanium detectors at high counting rates 
and higher energies [1]. Jordanov et al. developed the recur- 
sive algorithm of trapezoidal pulse shaping by using a con- 
volution method [2, 3]. The algorithm based on convolution 
method made trapezoidal pulse shaping possible by digital 
means. Stein et al. used moving window deconvolution tech- 
nique to realize trapezoidal and triangular shaping in digital 
signal processor (DSP) [4]. In addition, a real time digital 
pulse shaper, which was synthesized by a concave and a con- 
vex pulse shapes, was also implemented in programmable 
logic [5]. Also, z transform method was described to ob- 
tain the recursive algorithm of trapezoidal pulse shaping by 
Cosimo Imperiale [6]. Some simulations of the algorithm 
based on z transform method were discussed [7-10]. Re- 
cently, Trapezoidal shaper was employed in digital spectrom- 
eter for ballistic deficit correction, neutron-gamma discrimi- 
nation and pile up correction [11, 12]. Furthermore, typical 
shapers, including triangular, trapezoidal and cusp-like ones, 
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were generated in a new adaptive digital shaper which en- 
abled automatic adjustment of coefficients for shaping an in- 
put signal [13]. However, when the input signal is associated 
with noise, the baseline drift of trapezoidal pulse appears. It 
is especially obvious in the process of consecutive pulses. In 
this paper, we propose two methods to remove the baseline 
drift. The feasibility and accuracy of the methods are verified 
by simulations and experiments. The results show that the 
methods can remove the baseline drift in trapezoidal pulse 
shaping. 


Il. TRAPEZOIDAL PULSE SHAPING 


The trapezoidal pulse shaping algorithm was raised by z 
transform method [6]. As shown in Fig. 1, a trapezoidal pulse 
can be directly synthesized by Eq. (1): 


4 
v =X vi(t), (1) 
w=1 


where, v1(t) = (Vinax/ta):tu(t), vo(t) = —v1(t—t,)u(t—ta), 
u3(t) = —v1(t—t))u(t—t,), and v4(t) = vı (t—t.)u(t—te); 
ta is the rise time of trapezoidal pulse, tp — ta is the duration 
of flat top, te is the total width of pulse, and Vmax is the height. 

Equation (1) can be described as Eq. (2) according to z 
transform 


Vinax2(1 — 27" — 27M 4+ 27") 

Volz) = : 2 

where, na = ta/Ty, ne = ty/Ty, Me = te/T;, and T, is the 
sampling time of ADC. The input signal is defined as 

v(t) = Ae~/Tu(t) t >0, (3) 

where A is the height and 7 is the time constant. The function 

of the input signal in z transform notation is 
1 

1—dz-?’ 


Viz) =A (4) 
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where d = e~7s/7 


represented as 


. The transfer function of Fig. 1 can be 


owe) a- dz“")\(l—27-™)(1 — 27™)z71 
Ha Vi(z) ma(l — z1)? 
(5) 


(1 —dz-1)(1 — 2-™)(1- z7) t 


Vo(2)(1-271)? = (2): 


Na 
(6) 
Applying inverse z transform, we have the time domain 
output of vo: 


Nna — 1] — vi[n — m — 1] + viin — ne — 1] —d- [viln — 2] 


vijn — na — 2] — vi[n — ny — 2] + viln — ne — 2]]} - — - 
vojn] = viln] = 0 n<0 


A yin] = 2y[n — 1] — yn — 2] + Xn], (8) 
E 
2 where 
2 
= 1 
g X[n] = {z[n — 1] — z[n — na — 1] — z[n — m — 1] 
2 Na 
s 
g + z[n — ne — 1] — d - [a[n — 2] — x[n — na — 2] 
—al[n — w% — 2] + z[n — ne — 2]]}. 
0 1024 2048 3072 4096 
Time (ns) The recursive equations of Eq. (8) are as follows: 
Fig. 2. (Color online) Trapezoidal pulse shaping with baseline drift. yin] — y[n — 1] = y[n — 1] — y[n — 2] + X [n] 
yin- 1] -yin — 2] = yin — 2] - yin - 3] + X[n — 1] 
—2|— — 3| = =3 = —4|+Xjħn- 2 
The trapezoidal pulse shaping is implemented by Eq. (7). ylr — 2] — ln uns =y =a] =a 
It is suitable for processing ideal input signal without noise. ; 
Generally, the input signal is superimposed by Gaussian l 
white noise. An input signal with SNR = 30dB is sim- y[3] — yl2] = y[2] — yf] + XB) 
ulated by Eq. (3) with A = 2000 and r = 100. The re- y[2] — y[1] = y[1] — y[0] + X [2]. 
sult of trapezoidal pulse shaping with na = 150, np = 300, 
ne = 450, T; = 1 is shown in Fig. 2. The baseline drift of An equation only about y[n], yin — 1], y[1], y[0] and X [n] 
trapezoidal pulse appears when the input signal is associated is acquired by accumulating the equations above 
with Gaussian white noise. The simulations results indicate n 
that the lower SNR, the more serious baseline drift will be. : 
yl] — yin- 1] = yi1] — yl] + D> XT. 
i=2 
Ill. SIMULATIONS AND EXPERIMENTAL TESTS The general formula of y[n] is 
3 
A. Convergence conditions 2 
yin] — yit] = (n — 1)(y[1] — yl0]) + x12] + $ Xi] 
i=2 
Figure 2 shows that the algorithm based on z transform n—1 m (9) 
method is not convergent. To remove the baseline drift, a deepak 5 Xļ|i] + 5 Xi]. 
new method based on convergence analysis is used. Letting pa par 
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Output 2 
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Time (ns) 
Fig. 3. (Color online) Baseline drift removal based on the new 
method. 


Assuming that y[0] = y[1] = 0, one can simplify Eq. (9) 
into Eq. (10): 


yin] xay x+ ++ X. (10) 


vijn] is an equal ratio progression, with an equal ratio of 
q = e-\/7. It can be proved that X[n] is also an equal ra- 
tio progression and the equal ratio is g. Equation (10) is ex- 
pressed as 


yin] = X[2] + xBu- sia Ba- ) 
XB- 
l1—q : 
(11) 


Let y[n] = X[2] + Y, where 
Y=C-(-¢@)+---4+C-(1-—q™*)+C-(1-q), 
where C = X[2]/(1 — q), and then 


Y =Chn-2- P0- 7)/0-a)] aD 
q < last > 1. So, y[n] can be convergent when Y = 
0. It can be seen from Eqs. (8) and (12) that Y can be set 
to 0 by setting x|1] to 0. Figure 3 illustrates the result of 
trapezoidal pulse shaping with v;[1] = 0, with the same input 
signal as Fig. 2. Output 1 is used as a contrast. Output 2 is the 
trapezoidal pulse with the first input data being 0. Therefore, 
the new method can remove the baseline drift resulting from 
trapezoidal pulse shaping. 


Nucl. Sci. Tech. 26, 050402 (2015) 


B. Noise suppression 


The baseline drift of trapezoidal pulse is caused by noise 
accumulation. Filtering the original signal before shaping is 
another way of removing the baseline drift. Digital S-K filter 
performs well in signal processing with amplitude filtering 
and frequency filtering factors [14]. The true height of filtered 
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Fig. 4. (Color online) Digital S-K filter in trapezoidal pulse shaping. 


signal can be obtained by adjusting amplitude of the filtering 
factor properly. The algorithm of digital S-K filter is written 
as 


(k-(3—a)+2k?)-y[n—1]—k?-y[n—2]+a-2[n] 
1+k-(3—a)+k? 
n<0O 


. 8 
x[n] = 0 os 


m 
SES 
ai 
Io l 


where x[n] is the discrete input signal, y[n] is the output sig- 
nal, k is the frequency filtering factor and a is the amplitude 
filtering factor. The original signal in Fig. 2 is filtered by the 
digital S-K filter, and output of the filter is processed by call- 
ing Eq. (7) recursively. The baseline drift of trapezoidal pulse 
is removed with k = 5, a = 1.15, as shown in Fig. 4. 


Differential operation can also attenuate the noise. Equa- 
tion (6) can be rewritten as 


(1 —dz7!)(1— 27-™)(1 — z7™)z71 
na(1 — 2-1)? 


Vo(z) = Vi(z): . (14) 


Assuming V(z) = V(z)/(1 — 1/z)?, an improved recursive 
algorithm of trapezoidal pulse shaping is obtained: 


vin] = 2u[n — 1] — v[n — 2] + viln] 

Vo[n] = A {v[n — 1] — v[n — na — 1] — v[n — m — 1] + o[n — ne — 1] (15) 
—d-(v[n — 2] — v[n — na — 2] — v[n — m — 2] + vin — ne — 2])} 

voln] = v[n] = 0 n<0 
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Fig. 5. (Color online) Trapezoidal pulse shaping based on the im- 
proved recursive algorithm. 


Equation (15) is the improved recursive algorithm of trape- 
zoidal pulse shaping without baseline drift [10]. Trapezoidal 
pulse can be implemented by calling Eq. (15). Figure 5 shows 
that the improved algorithm can remove the baseline drift in 
Fig. 2. 


C. Measured pulses shaping 


In order to verify feasibility of the methods, pulses tests 
were carried out. The time constant of pulses is 3.2 us, de- 
tected by silicon drift diode (SDD) detector and digitized by 
ADC at 20 MSPS, with T; = 50 ns and 7 = 3.2 us according 
to Eq. (7). A trapezoidal pulse with ta = 4 us, t = 4 us and 
te = 8us corresponding to na = 80, np = 80 and n, = 160 
is implemented by Eqs. (7) and (15). The experimental tests 
of different methods are given in Fig. 6. 

The results indicate that both the new method and the im- 
proved recursive algorithm can remove the baseline drift in 
trapezoidal pulse shaping. For the purpose of signal filtering, 
digital S-K filter can also be used in trapezoidal pulse shap- 
ing to remove baseline drift. However, it is complex to realize 
digital S-K filter in FPGA. Besides, it is also difficult to se- 
lect the values of k and a. Therefore, the new method and 
the improved recursive algorithm are recommended in prac- 
tical applications. The correlation between the new method 
and the improved recursive algorithm is shown in Fig. 7. It 
indicates that the improved method is suitable for trapezoidal 
pulse shaping without baseline drift. 


IV. CONCLUSION 


In this paper, the recursive algorithm of trapezoidal pulse 
shaping is derived by z transform method and the approaches 
to baseline drift removal are discussed. Setting the first in- 
put data to O and digital S-K filter are used to remove the 
baseline drift due to the accumulation of noise. The results 
of experimental tests demonstrate that the methods can re- 
move the baseline drift existing in trapezoidal pulse shaping. 
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Fig. 6. (Color online) Pulses shaping with (a) the first input data be- 
ing 0, (b).digital S-K filter processing and (c) the improved recursive 
algorithm. 
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Fig. 7. (Color online) Correlation between the new method and the 
improved recursive algorithm. 


050402-4 


NEW METHODS TO REMOVE BASELINE... 


A comparison between the new method and the improved re- 
cursive algorithm is also carried out. It shows that the new 
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